Low-Energy Quasiparticles in Cuprate Superconductors: A Quantitative Analysis 
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A residual linear term is observed in the thermal conductivity of optimally-doped Bi2Sr2CaCu2 0s 
at very low temperatures whose magnitude is in excellent agreement with the value expected from 
Fermi-liquid theory and the d-wave energy spectrum measured by photoemission spectroscopy, with 
no adjustable parameters. This solid basis allows us to make a quantitative analysis of thermo- 
dynamic properties at low temperature and establish that thermally-excited quasiparticles are a 
significant, perhaps even the dominant mechanism in suppressing the superfluid density in cuprate 
superconductors Bi2Sr2CaCu2 08 and YBa2Cu307. 

PACS numbers: 74.70.Tx, 74.25.Fy 
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The superconducting order parameter of the archety- 
pal high Tc compounds YBa2Cu307 and Bi2Sr2CaCu208 
is widely agreed to have d-wave symmetry, yet there is no 
consensus on the correct theoretical description of their 
superconducting state properties, let alone those of the 
metallic state. A fundamental issue in the current de- 
bate is the nature of the electronic excitations in tjiese 
systems, and vifhether long-lived quasiparticles existtJ or 
not.a Another debate concerns the dominant mechanism 
responsible for the thermal suppression of the superEuid 
density, whether it, be d-wave nodal quasiparticlesQ or 
phase fiuctuationSjQ for example. One way to shed new 
light on these issues is to go beyond the usual qualita- 
tive temperature dependence of physical properties, and 
look closely at their magnitude. Our specific approach 
is to examine quantitatively the basic thermodynamic 
and transport properties of these two superconductors 
within a Fermi-liquid description of d-wave quasiparti- 
cles, grounded in a spectroscopic measurement of the en- 
ergy spectrum, and see whether a consistent description 
at low energies can be achieved. 

The d^2_y2 gap function goes to zero at four nodes 
along the k^ = iky directions, producing a conelike 
quasiparticle excitation spectrum at low energies: 
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where vf and V2 are the energy dispersions, or quasi- 
particle velocities, along directions normal (|j ki) and 
tangential (|| k2) to the Fermi surface, at each node. 
This spectrum is associated with the two-dimensional 
Cu02 plane that is the fundamental building block of 
all cuprates. It neglects any possible dispersion in the 



third direction (along the c-axis), as well as excitations 
associated with the one-dimensional CuO chains found in 
some crystal structures, notably in YBa2Cu307_i5 (along 
the 6-axis). 

This simple spectrum gives rise to a quasiparticle den- 
sity of states which is linear in energy: 
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which in turn leads to a T^ dependence of the elec- 
tronic specific heat and a linear T dependence of the su- 
perfluid density, for example. In a realistic treatment, 
one needs to include the effect of impurity scattering and 
electron-electron interactions. One usually accounts for 
the former in terms of a single, isotropic scattering rate, 
parameterized by an impurity bandwidth 7. At ener- 
gies below 7, known as the "dirty" limit, one expects a 
profound modification of the density of states which ac- 
quires a residual finite value A^(0) ex 7. At energies well 
above 7, in the "clean" limit, N{E) ex E and one recovers 
many of the straightforward temperature dependences. 
Going beyond this, Durst and Lee recently included ver- 
tex corrections, which arise because of the anisotropy 
of scattering in a d-wave superconductor. Q The impor- 
tance of Fermi-liquid corrections has also been empha- 
sized, whereby electron-elefjlron interactions renormalize 
the normal fluid density.EfEli3 

In this paper, we use a measurement of the in-plane 
thermal conductivity at very low temperature to ex- 
tract a value for the ratio vp/v2 in optimally-doped 
YBa2Cu307-d- (YBCO) and Bi2Sr2CaCu208 (BSCCO). 
We proceed to show first that for BSCCO this ratio is in 
excellent agreement with the values of vp and V2 mea- 



sured separately by angle-resolved photoemission spec- 
troscopy (ARPES). We then use the ratio to calculate 
the drop in superfluid density within a Fermi-liquid de- 
scription and compare this with the experimental results 
obtained from penetration depth measurements. Finally, 
we extend our quantitative analysis to include specific 
heat measurements in YBCO. From the overall analy- 
sis, we conclude that the superconducting state of the 
cuprates is well described by Fermi-liquid theory, at least 
at low energy and optimal doping. 

Thermal conductivity The thermal conductivity of 
YBa2Cu306.9 and Bi2Sr2CaCu208 was measured using a 
steady-state method, described elsewhere.0 The samples 
were single crystals grown via standard flux techniques 
and oxygenated so as to obtain the maximum Tc (opti- 
mal doping), quoted in Table I. Both crystal structures 
are made of Cu02 planes stacked along the c-axis, with 
a density 30% higher in YBCO, due to its lower average 
interplane spacing, given in Table I. The conductivity 
of YBCO was measured along the a-axis in untwinned 
crystals, so as to avoid the contribution of CuO chains. 

In Fig. 1 we present the low-temperature thermal con- 
ductivity, K, of BSCCO, c«d compare it with that of 
YBCO obtained previouslyJa By plotting n/T vs T^, we 
can separate the linear quasiparticle term from the cu- 
bic phonon term (see Ref. R). A finite residual linear 
term k^/T (the value of k/T as T ^ 0) is observed, of 
similar magnitude for the two cuprates, given in Table I. 
(Note that tke value for YBCO is an average over sev- 
eral samples,! only one of which is displayed in Fig. 1.) 
The error bar on these numbers is approximately ±20%, 
arising about equally from the uncertainty in the extrap- 
olation and in the geometric factor of each sample. 

Calculations for the transport (QOieat by d-wave quasi- 
particles in two dimensions givefl'u 
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where n/d is the stacking density of Cu02 planes. The 
residual conduction is due to a fluid of zero-energy quasi- 
particles induced by the pair-breaking effect of impurity 
scattering near the nodes in the gap, and it is indepen- 
dent of impurity concentration. This universal character 
of kq/T waSj-demonstrated explicitly for both YBCOQ 
and BSCCO. y 

Durst and Lee recently showed Eq. (3) to be valid even 
when vettex and Fermi-liquid corrections are taken into 
account p so that unlike charge transport, heat transport 
is not renormalized by either correction. The universal 
character and the absence of renormalization make ther- 
mal conductivity a privileged probe of the quasiparticle 
spectrum in a d-wave superconductor, providing a simple 
and direct measurement of vf/v2 in the cuprates. From 
the measured kq/T and the known values of n/d (see 
Table I), we obtain: 
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FIG. 1. Thermal conductivity divided by temperature vs 
T of YBa2Cu307 (squares) and Bi2Sr2CaCu208 (circles), at 
optimum doping. The lines are linear fits to the data below 
130 mK, with extrapolated values given in Table I. 
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with an uncertainty of about ±20%. 
ratio for YBCO is twice the value of '-- 
the literature.) 
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ARPES Angle-resolved photoemission spectroscopy 
has established the existence of a Fermi surface in 
YBCOy and BSCCOBfl and revealed directly the 
k-dependep£e_ of the gap characteristic of dx^-y^ 
symmetry£3'E3 In BSCCO, the nodes are along the (0,0) 

to (7r,7r) direction, at k = kp = 0.74 A Ji3 where the 
energy-has a dispersion along ki given by vp — 2.5 x 10^ 
cm/slj As for the dispersion along k2 (or 0), Mesot 
et alx3 were recently able to extract S = \dA/d(f>\node, 
the slope of the gap at the node. For a crystal near 
optimal doping {Tc — 87 K), they obtain S — 60 ± 5 
meV (= 1.7Amax, where A,„ax is the gap maximum, at 
cf) = 0). This yields V2 = S/hkp = 1.2 x 10^ cm/s, so 
that Vp/v2 = 20. 

A hotly debated question is whether the excitations 
in the vicinity of the Fermi surface, in particular along 
the diagonals, |-Gan be treated as the usual Landau/BCS 
quasiparticles.ui3 The excellent quantitative agreement 
we find between the spectroscopic and the transport 
measurement of vp/v2 in BSCCO strongly validates a 
Fermi-liquid description of the superconducting state in 
cuprates, at least at low energies. 

Unfortunately, ARPES measurements in YBCO have 
been less successful so far. The Fermi surface is mO|m 
complicated, with bilayer splitting of the plane bandsEiJ 



TABLE I. Comparison of YBa2Cu306.9 and Bi2Sr2CaCu208 at optimal doping. In YBCO, all directional properties are 
for the a-axis (no chain contribution), d/n is the average separation between Cu02 planes-stacked along the c-axis. The 
zero-temperature penetration depth A(0) was measured by far-infrared reflectivity in YBC0E3 and bW^ magnetisation in 
BSCCO.Ea The Fermi velocity vp and momentum kp were obtained from angle- resolved photoemission.Eil'll3 vpjv^ is the ratio 
of quasiparticle velocities, obtained via Eq. (3), using the residual linear term kq/T measured in the thermal conductivity at 
T — > 0. S = dl\{ij>) / d(j) = hkpV2 is the slope-pM^he gap at the node calculated using kp, vp and vpjv-i. Amax is the gap 
maximum as seen in c-axis tunneling by STM,cjc3 with S — /iAmax. The linear drop in superfluid diessity with temperature is 
expressed as A^(0)dA~^/dr, obtained from the penetration depth measured at microwave frequencies .Ej'Ell c? is the Fermi-liquid 
correction computed from Eq. (7), using the measured values of d\~^ jdT and kq/T. 
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and an added band for the CuO chains. When aver- 
aged over the two plane bands, the band crossing and 
dispersion at the Fermi energy along the (0,0) to (tt, tt) 
direction are close to those quoted above for BSCCO, 
namely vp — 2.5 x 10^ cm/s and kp — 0.8 A~^,Eil albeit 
with greater uncertainty. The gap structure has not yet 
been resolved with sufficient resolution to provide a mea- 
surement of V2. The thermal conductivity data may be 
used instead: with vp — 2.5 x 10^ cm/s, Eq. (5) yields 
V2 — 1.8 X 10^ cm/s. This implies that the slope of 
the gap at the node in YBCO is 1.5 times larger than in 
BSCCO, with 5' — hkpV2 — 95 meV, in contrast with ev- 
idence from STM measurements of c-axis tunneling that 
the gap maximum in YBiC.Ci.is smaller thaip_in BSCCO, 
namely A^ax ^ 20 meVBEj vs ~ 40 mevN. This sug- 
gests a strikingly different angular dependence of the gap 
function, with a ratio of slope to gap maximum 3 times 
larger in YBCO (see Table I), under the assumption that 
Vp is the same in both materials. 

Superfluid density As the temperature is increased 
from T = 0, the thermal excitation of nodal quasipar- 
ticles causes the normal fluid density, p„(T), to grow 
linearly witfut^^ijipjerature. In the clean limit, at low 
temperaturem 
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where m is the mass of the carriers jEiid a^ is the Fermi- 
liquid correction for charge currents £3 

A linear temperature dependence of p„ (T) is a charac- 
teristic feature of most cuprate superconductors, as re- 
vealed through measurements of the penetration depth 
A(T), via the relation ps{T)/m — ps{0)/m — pn(T)/m = 
c^/47re^A^(T). From the data of-Jferdy and co-workers 
on untwinned crystals of YBCOt2[E3 - again taking the 



a-axis results to avoid chain contributions - and from 
measurements by Waldram and co-workers in BSCCO, Eil 
one finds the slope of A^(0)/A^(T) at low temperatures, 
given in Table I for optimal doping. Combining Eqs. (3) 
and (6), we can then solve for the Fermi-liquid correction, 
via: 



dA-^(r) 

dT 



^2.93 X 1013 ^ a' 



(7) 



with A in meters and kq/T in W K ^ m ^. Using the 



values for A(0) quoted in Table I, we get: 
BSCCO : a^ = 0.43 



YBCO 



= 0.46 



(8) 
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In other words, the observed drop in superfluid den- 
sity is about 2 times weaker than expected from a cal- 
culation neglecting interactions, as shown graphically for 
YBCO in Fig. 2. Since electrons in cuprates are highly 
correlated, a renormalization by a factor of 2 seems en- 
tirely plausible. The fact that it is comparable in the 
two compounds is not unexpected, given that the Fermi 
velocities, themselves renormalized by interactions ,13 are 
comparable. We stress that an estimate of the Fermi- 
liquid correction to Ps(T) does not require a separate 
knowledge of vp and V2, and heat conduction - unlike 
heat capacity, for example - can provide directly the ap- 
propriate combination of the two parameters, i.e. th 
ratio. Note that in the case of BSCCO, Mesot et all 
were the first to report an estimate of the renormaliza- 
tion factor, based on their ARPES data. (The fact that 
they obtain a slightly different value, namely a^ = 0.32, 
is due to their use of different penetration depth data, 
which we consider to be less reliable because restricted 
to temperatures above 17 K.) 
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FIG. 2. Temperature dependence of the superfluid density, 
normalized to unity at T = 0, for optimally-doped YBCO 
(a-axis) (from Ref. KOf). The lines are the expected low tem- 
perature behaviour calculated from Eq. (7) with (a^ 7^ 1) and 
without (a^ = 1) Fermi-liquid interactions. 



Given the numbers that emerge from the analysis, it 
seems fair to conclude that the thermal excitation of 
quasiparticles is a significant, perhaps even the domi- 
nant mechanism in suppressing the superfluid density 
of these two cuprate superconductors. It is interesting 
that electron-electron interactions appear to be such as 
to weaken this process. 

It is perhaps worth noticing that although the density 
of superfluid at T = (ex A~^(0)) is 1.7 times higher 
in YBCO, the normal fluid density grows at exactly the 
same rate in both compounds at low temperature. This 
would naively suggest that Tc should be much higher in 
YBCO, while it is in fact not very different (5% higher). 
This is because Ps{T) acquires a much stronger down- 
ward curvature near Tc in YBCO. Therefore, a major 
difference must develop at higher temperatures. Part of 
the answer must come from the very different curvature 
of the gap function away from the nodes, at high ener- 
gies, as parametrized by /i = S'/Amax being much larger 
in YBCO (see Table I). 

Specific heat We complete our quantitative analysis 
by looking at the electronic specific heat, which is simply 
derived from Eq. (2), in the clean limit: 
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where C(3) ~ 1.20. Extracting this electronic contribu- 
tion from the total specific heat has been a controversial 
exercise, since the data can be fitted equally well with- 
out a r^ term. In YBCO at optimal doping, the value 



quoted in the literature is 0.1 mJ K~'^ mole~^,cJS albeit 
with a ± 60% uncertainty. Using vp = 2.5 x 10^ cm/s and 
V2 — 1.8 X 10^ cm/s, Eq. (10) gives 0.065 mJ K^-^ mole~^ 
- a value within the experimental uncertainty. 

An alternative approach is to extract V2 from the field 
dependence of the specific heat. The Doppler shift of 
quasiparticle states near the nodes in the presence of the 
superfluid flow around vortices leadS|-tp an increase in 
the specific heat proportional to \/h!E3 In terms of the 
nodal spectrum, the magnitude of the effect is related 
only to the slope of the gap at the node, V2- In the 
clean limit, the electronic specific heat of Cu02 planes 
(per unit volume), calculated by averaging the effect of 
the Doppler shift over a single vortex-lattice unit cell, is 
given by nil 
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where $0 is the flux quantum, and a is a vortex-lattice 
parameter of order unity. Such a yfj dependence has 
been seen in measurements on YBCCO'E3 and the co- 
efficient, of magnitude 0.9 mJ K~^ mole^^ T^^/^, de- 
termined with greater accuracy (±10%) than the corre- 
sponding zero-field T^ term. Using V2 = 1.8 x 10^ cm/s 
and a = 1 in Eq. (11) gives 0.6 mJ K~^ mole^^ T^^/^. 

It is clear that both aspects of the specific heat data for 
YBCO are in reasonable quantitative agreement with our 
thermal conductivity result for vp /v2 and the value of vp 
from ARPES. In a refined treatment, one would take into 
account the contribution of CuO chains to the density of 
states (and hence to the specific heat). In this respect, it 
is interesting that Junod and co-workers extract T^ and 
vi? coefficients for sio-overdoped crystal of YBCO which 
are somewhat larger:Ea 0.20±0.05 mJ K~'^ mole~^ and 
1.3±0.1 mJ K~^ mole~^ T^^/^. It is not unreasonable to 
attribute this increase to a larger chain density of states, 
such as would explain the decreasing a— h anisotropy in 
the linear temperature drop of the superfiuid density ob- 
served in crystals of YBCO as one moves from overdoped 
to underdoped.c3 

In conclusion, we have provided a quantitative anal- 
ysis of low-temperature data for the cuprate supercon- 
ductors YBCO and BSCCO at optimal doping which 
compared results from our thermal conductivity mea- 
surements with existing results from ARPES, microwave 
penetration depth and specific heat. Within a Fermi- 
liquid model of d-wave quasiparticle excitations with 
interactions, we find all data consistent with a single 
set of parameters. The Fermi velocity and the Fermi- 
liquid renormalization of charge currents are found to be 
roughly the same in both compounds (as is Tc), with 
Vp ~ 2.5 X 10^ cm/s and a^ ~ 0.4 — 0.5, whereas the 
slope of the gap at the node is about 1.5 times steeper in 
YBCO. 

In particular, the thermal excitation of quasiparticles 
emerges as a sufficient mechanism for suppressing the su- 
perfiuid, and there is no clear evidence for a significant 



contribution from phase fluctuations at low temperature, 
at least at optimal doping. The success of a Fermi-liquid 
description for the low-temperature properties should not 
be taken to mean that the normal state of cuprates is a 
Fermi liquid. Nor should it be viewed as supporting a 
BCS theory of the superconducting state, given the fact 
that, for example, the thermal excitation of quasipartk 
cles does not have the expected impact on the gap itselLS 
which remains undiminished up to high temperatures .Ii3 
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